UKQCD Collaboration and University of Edinburgh 



LATTICE QCD WITH SUPPRESSED HIGH MOMENTUM 
MODES OF THE DIRAC OPERATOR 

ARTAN BORigi 



Abstract. I define lattice fermions in five Euclidean dimensions and the cor- 
responding effective theory in four dimensions. The main properties of these 
theories include the suppression of high momentum modes of the lattice Dirac 
operator and their ability to continuously interpolate between quenched and 
dynamical fermions. In particular, the standard formulation of lattice QCD 
can be viewed as a limiting case of the theory. 
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Lattice simulations are an indispensable tool in understanding the strong force. 
Since its formulation by Wilson, 1974 1, lattice Quantum Chromodynamics (QCD) 
has progressed into a separate discipline. 

Nevertheless, simulations of lattice QCD are still away form the precision tests 
that one would like to see. Clearly, faster computers and simulation algorithms are 
needed. 

Recently, the issue of 'ultraviolet (UV) slowing down' has attracted particular in- 
terest [Irving et al, 199S, Duncan, Eichtcn and Thackcr, 1999, de Forcrand, 1999, 
Peardon, 2001, A. Hasenfratz fc Kncchtli, 2001| . These studies try to address al- 



gortmically large fluctuations of the high end modes of the fermion determinant. 
The goal of these algorithms is to separate the UV modes and to focus on the 
silumation of the physically interesting infrared modes. In this paper I show that 
all this computational effort can be saved if one improves UV properties of the 
lattice Dirac operator in the first place. 

Much of U V slowing down is thought to c ome from the non-smoothness of the 
gauge fields | A. Hasenfratz fc Kncchtli, 2001 . The effects of non-smooth gauge 



fields are mostly observed in the large eigenvalues of the fermion determinant 
Irving et al, 1998 1 . Recently, | Duncan, Eichten and Thacker, 1999 1 proposed a strat- 
egy to accelerate fermion simulations by computing directly the 'infrared' eigenval- 
ues and attaching the 'ultraviolet' ones by approximate actions and the multiboson 
method, de Forcrand, 1999 proposed an algorithm which 'filters' these fermion 
modes. Another strategy is the inclusion of ultraviolet modes in a multibosonic 
fashion which results in faster algorithms [Peardon, 2001]. A direct smoothing ap- 
proach is also possible by smearing techniques [A. Hasenfratz & Kncchtli, 2001]. 

In spite of the recent progress there is no unifying view how to deal with the UV 
slowing down. In this paper I propose an improved formulation of lattice fermions 
which 'quenches' fermion eigenvalues beyond a given eigenvalue level. The paper 
proves rigorously the existence of such a field theory. 

In the following section I discuss the need to deal with the problem of non-smooth 
lattice gauge fields. In section 2 I give a model of five dimensional fermions with 
suppressed high fermion modes. I then define the theory in four dimensions and 
prove its the field theoretic properties in section 3. Finally, in section 4 I draw the 
conclusions. 



1. Difficulties with lattice fermions 

1.1. Basic defin itions. The lattice regularization of gauge theories was defined 
by jWilson, 1974 1. 

A fermion field on a regular Euclidean lattice A is a Grassmann valued vector 
ip(x), x — {x^,)! = 1, ... ,4} G A which carries spin and color indices. The first 
and second order differences are defined by the following expressions: 



d^{x) = j^[ip(x + ae^j - ip(x - ae^)} 

dfa(x) = -^[i>{x + ae M ) + tp(x - ae M ) - 2i/j(x)] 
where a and e„ are the lattice spacing and the unit lattice vector along the coor- 



dinate [i — 1,... ,4. Let U(x)^ G 



-3x3 



be an element of the SU{i) group, the 



oriented link connecting lattice sites x and x + ae^. 
defined by: 



Then covariant differences are 
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V^(x) = ^[U(x) fl ip(x + ae^) - W(x - ae^)^ip{x - oe M )] 

A M -0(x) = -^[U(x) IJ ,ip(x + ae M ) + W(x- ae^^x - ae M ) - 2ip(x)] 
The Wilson-Dirac operator is a matrix operator Dw G C NxN defined by: 

4 

(1-1) £> w =ml + ^( 7M V At -|A M ) 

where 1 is the identity matrix and m the bare mass of the fermion; {7^ G C 4x4 , /i = 
1, . . . , 4} is the set of anti-commuting and Hermitian gamma-matrices of the Dirac- 
Clifford algebra. The fermion lattice action is defined by: 

(1.2) S f = $(x)Dw(x,yWy) 

whereas the gauge action is given by: 

(1.3) S 9 = 4^(1-1%) 

9 p 

The sum in the right hand side is over all plaquettes V on the lattice. Up is 1 x 1 
Wilson loop and g is the bare coupling constant of the theory. 

The basic computational task in lattice QCD is the evaluation of the partition 
function given by: 

(1.4) Z QCD = f a H {U)a{>>p,i>)e- s f- s ° 



where <jh(U) and a(ip,tp) denote the Haar and Grassmann measures respectively. 
The computing problem has 0(e N ) complexity (i.e. it is NP-hard) and one has to 



resort to stochastic estimations of the right hand side (1.4). In fact, integration 



over the Grassmann fields can be performed exactly to give: 

Zqcd = J <Jh (U) det %e _s » 

1.2. A measure of non-smoothness of lattice gauge fields. I will use a linear 
function of plaquette to characterize the degree of the non-smoothness on the lattice. 
Let P{U) be the average lxl Wilson loop over all lattice plaquettes. Then I define 
the following real and positive function: 

(1.5) $(Z7) = 1 - p(U) = 1 - (U P )p 

Smooth gauge fields are characterized by small values of $>(U). Typical values 
from lattice simulations arc §>(U) € (0.4,0.5), which is a clear indication of non- 
smoothness of the lattice gauge fields. 

To see the influence of the non-smoothness on a typical gauge invariant lattice 
operator, I consider eigenvalue perturbations of the lattice Dirac operator. Let 
X(U) be an eigenvalue of the matrix A{U) = D^Dy/- A classical result from the 
eigenvalue perturbation theory states that variation of an eigenvalue 8\{U) under 



matrix perturbation 5A(U) is bounded by [Golub fc Van Loan, 1989|: 



(1.6) \5X(U)\ < \\SA(U)\\ p , peN 
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The following result relates the 2-norm of the matrix A(U) to the degree of the 
non-smoothness $(f ) on the lattice (fbq): 



(1.7) \\8A(U)\\ 2 < ci + cay/tyJJ) 

where c\ and C2 are positive constants. This result is proven in Appendix A. 
Together with the bound (|l.6[) it gives: 



(1.8) \6\{U)\ < Cl +c 2 yMu) 



The above inequality (1.8) states that the eigenvalue fluctuations on a lattice 
with smooth gauge fields are likely to be smaller than those on a lattice with non- 
smooth gauge fields. 

However, as present computing resources do not allow a close approach to con- 
tinuum limit, it is important to look for formulations and algorithms with reduced 
effects of lattice non-smoothness. 

Numerical results of |Duncan, Eichten and Thacker, 1999 1 indicate that large 



eigenvalues of the lattice Dirac operator are merely lattice artifacts. It is a well- 
known fact that cutoff modes are poorly represented on the lattice. But one may 
not simply exclude them from simulations since the existence of the theory may be 
compromised. The strategy followed in this paper is the suppression of the high 
fermion modes of the lattice theory. As shown below it is possible to model a 
fermion theory with the reduced appearance of these modes. 

2. Modeling cutoff modes: Wilson fermions in 4+1 Euclidean 

dimensions 

Recent progress with chiral fermions on the lattice has shown that a theory of five 



dimensional fermions can be a useful modeling tool (for a review see [Kikukawa, 2001 
Domain wall boundary conditions along the fifth dimension provide a kinematical 
model for QCD with chiral fermions which are "localized" on the surface of a five 
dimensional world. The theory in five dimensions can be viewed as a fermion sys- 
tem propagating along the fifth Euclidean dimension with its dynamics generated 
by a certain Hamiltonian operator Ti. 

Let c*(x) and c(x) be creation and destruction operators in the Fock space which 
satisfy the anti-commutation relations: 

[J(x),c(x)]+ = 1 

They carry spin and color indices, which are not explicitly shown for clarity. c{x) 
acts on the bare vacuum state which they annihilate. 

I define the Hamiltonian operator of the fermion system by the bilinear form: 

"H = ^ J(x)H w (x,y)c(y) 

where H\y is the Hermitian lattice Dirac operator in four dimensions given by: 

Hw = IsDw 

Let L5 be the lattice size in the fifth dimension or the "the inverse temperature" 
of the quantum statistical system with the partition function Z, given by: 

Z = Tr er L5U 
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To compute the trace of an operator in the Fock space I use the standard tech- 
nique of the Grassmann coherent states with antiperiodic boundary conditions in 
the fifth Euclidean coordinate. Since Tt is quadratic it is easy to show that: 

Z = det(l+ e - L&Hn ) 

I define the measure density of the non-trivial fermion theory on the lattice by 
the following equation: 

u>{H w ) = iTr Oe- L * H 

where O is a function of the creation and annihilation operators. I choose it such 
that the right hand side is given by: 

u)(H w ) = - det(l - e - LsHw ) 
Such an operator exists as it is shown in the framework of domain wall fermions 



Furman, Shamir, 1995 



Hence the resulting density can be written as: 

(2.1) w(H w ) = dct(tanh y H w ) 

This suggests that an effective theory in four dimensions may be defined by the 
following lattice Dirac operator: 

(2.2) D = ^ stanh ^K 

a \l 

where a is the lattice spacing of the four dimensional lattice and [i > is a dimcn- 
sionless parameter. It is clear that for small lattice spacing this operator approaches 
the Wilson Dirac operator and hence has the correct continuum limit. 

In order to give to the Tr O operator a precise meaning, I define a fermion theory 
in five dimensions by the following action: 

(2.3) = $&D$1>W 

Here the five dimensional fermion field ip^ satisfies periodic boundary conditions 

(5) 

in all directions and D w is defined by: 

D*§ =D W + 75 V 5 - y A 5 



It can be shown that: The measure density of the five dimensio nal theory (2.5) 



is proportional to the measure of the effective theory defined by eq. (2.1). The proof 
is given in Appendix B. 

This result suggests that the four dimensional lattice theory with Wilson fermions 
can be approached by the 'high temperature' limit of a theory with Wilson fermions 
in five dimensions. Thus, it is natural to choose the length of the extra dimension 
to be proportional to the lattice spacing in four dimensions. Dimensional reduction 
is then realized by taking the continuum limit of the theory. Furthermore, the 
theory allows the introduction of a dimensionless parameter fj, which can be used 
to suppress the high momentum modes of the fermion theory to a prescribed level, 
i.e. \i can be viewed as a dimensionless 'temperature'. A 'cold' theory would then 
correspond to the quenched approximation, whereas a 'hot' one would be identical 
to the Wilson theory. 
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In the next section I will give the basic properties of the dimensionally reduced 
effective theory. 

3. A FERMION THEORY WITH SUPPRESSED CUTOFF MODES 

It is clear that the above theory defined in 4+1 dimensions has all desired prop- 
erties of a field theory: it is local, unitary and gauge invariant. 

This is not obvious for the dimensionally reduced theory with a lattice Dirac 



operator given by (2.2) 



u aHw 

3.1 £>=^-75tanh - 

a [i 

For free fermions the lattice Dirac operator has a momentum space representation 
given by: 

(3.2) D<p) = H 75 tarm^M 

a /i 

with p = {p vi v — 1,... ,4} being the four-momentum vector. The momentum 
space Wilson Dirac operator is given by: 

4 



a 

v=l 



whereas its square is given by: 

4 4 

a 2 b w {p)^I) w {p) = sin2 a Vv + ~ cos2 a P^ 2 

3.1. Locality. Whilst D\y links only nearest neighbour lattice points, D will be 
a full matrix. A full matrix can be considered essentially local if it is dominated 
by matrix elements which link lattice points that are close to each other. For 
example, this will be the case if the magnitude of decays exponentially with the 
distance \i — j\. This will be considered as a sufficient condition in the following for 



locality (see also [Hernandez et al, 1999 ) . Since D(p) is analytic and 27r-periodic 



for fi > 0, then its Fourier transform falls off exponentially at large distances (see 
Luscher, 199§| ] for a similar argument). Therefore, D is a local operator in the 



above sense. 

The locality of the fermion theory in a gauge field background is treated in 
Appendix C. In particular, it is shown that if the Wilson Dirac operator is singular 
the locality of the theory is guaranteed solely by the positivity of /i. 

3.2. Unitarity. For unitarity it is sufficient to show that the lattice operator leads 
to non-negative energy spectrum with non- negative norm of eigenmodes. To do this 
I define a positive function / in terms of the real variable z: 

f(z) — — tanhz, z ^ 

z 

Since tanhz is an odd function of z, one can easily show that the right hand side 
is in fact a function of z 2 only. Therefore, I can define a function / such that: 

, f(z 2 ) = itanhz, z ^ 

[6 - 6) /(0) = 1 
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This way, I may write: 

tanhz = zf(z 2 ), z e H 

Using the definition of an operator-valued function the lattice Dirac operator 
( |2.2| ) takes the form: 

(3.4) D = Dwf[{ ^K ) 2 ] 

M 

Now note that the matrix function /(.) is non-singular. Hence the poles in the 
fermion propagator are identical to those in the Wilson theory and the resulting 
theory is characterized by a real energy spectrum. Moreover, since /(.) is positive 
definite the norm of energy eigenmodes remains positive. 

3.3. Perturbation theory. The fermion propagator is given by the inverse of the 
expression (|]^). As usual, gauge fields are parametrized by su(3) elements: 

(3.5) U{x), t = e ia ^ A ^, A{x)„ £ su(3) 

and the Wilson operator is written as a sum of the free and interaction terms: 

Dw — D w + D w 

The splitting of the lattice Dirac operator is written in the same form: 

D = D a + D 1 , £>° = ^75tanh^ 
a \x 

where the interaction term has to be determined. This can be done by expanding 
D in terms of aj \x: 

(3.6) D = D w ^ + c 1 (^r+c 2 (^r + -..} 

where c\, C2, . . . are real expansion coefficients. 

Calculation of D 1 is outside of the scope of this paper. In fact it is an easy task 



if one stays with a finite number of terms in the right hand side of (3.6). Also, 
the number of terms can be minimized using a Chebyshev approximation for the 
hyperbolic tangent. [] 

3.4. Fixing /i. Tuning /i and then fixing it to a certain value is essential in spec- 
ifying the level of UV suppression and the fermion theory as such. On the prac- 
tical side one should know how to choose the value of /z such that only the in- 
frared modes are included. This can be done by computing explicitly the eigen- 
value density of the input operator Hyy and identify the threshold bet wee n the 



physical modes and the tail of the distribution. Follo wing definition (LI) and 



the technique of | Edwards, Heller and Narayanan, 1998 one can approximate the 



density of zero eigenvalues p(0, m) of —Hy/ at a given bare quark mass m. Fig. 
1 shows four example plots of low lying eigenvalue densities from the paper of 
(Edwards, Heller and Narayanan, 1998 1 . At small m the density is zero and then 



jumps at a non-zero value at the critical mass m cr where the Wilson operator be- 
comes singular. Then as m increases one can identify a threshold to m where the 
density approaches a plateau. It is this plateau where the UV effects dominate the 



I would like to thank Joachim Hein for discussions related to lattice perturbation theory. 
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Figure 1. Unnormalized density of zero eigenval- 
ues /?(0, to) of —H\y as a function of the bare quark 
mass to on 50 quenched configur ations as computed by 

Figure 



Edwards, Heller and Narayanan, 1998 1 . Figure taken from 
the above paper using its lanl.arXiv.org e-Print archive, 
hep-lat/9802016| version. 



eigenvalue spectrum. According to this heuristics one can suppress the eigenvalues 
beyond to^ and the value of /i can be determined by: 

(3.7) fi = rrifj, - m cr 
From Fig.l one can estimate fi « 0.3. 

3.5. Example of UV suppression. The effective action of a fermion theory can 
be written as: 

(3.8) S cS = trA 
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Figure 2 . Upper panel: Distribution of the noisy estimator X of 
the Wilson (left) and UV-sppressed (right) effective action. Lower 
panel: The same as in the upper panel but for the variation AX 
between two configurations. Lattice parameters: 12 3 x 24, (3 = 5.9, 
m = -0.869 (k = 0.1597), n = 0.2. 



where A = f(H^) is a Hermitian matrix and f(s) is a real and smooth function of 
s e M + . To see the effects of UV suppression one can compute the change of the 
effective action between two background SU (3) gauge fields. Since it is difficult to 
compute the trace directly one can use the noisy estimators of the type X = z T Az 
where z is a Z 2 niose vector. It is easy to show that the random variable X has 
expectation value E(X) = tvA and variance: 

VarpsT) =2^(Rea y ) 2 
Note that bilinear forms of the type z T Az can be computed using Lanczos based 



methods as described by |Bai et al, 1996, Cahill et al, 1999, Borigi, 2002b]. I have 
used the Lanczos algorithm as it appears in the paper of Borigi, 2002b |. Fig. 
2 shows the distribution of X and its variation AX between two configurations 
for Wilson and UV-suppressed fermions, i.e. f(s) — log y/s,s > and f(s) = 
logtanh v^s/Mj s > respectively. The figure shows clearly that changes AX in the 
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effective action estimator are reduced by an order of magnitude for UV-suppressed 
fermions as compared to Wilson fermions. 

4. Concluding remarks 

In this paper a lattice theory with suppressed cutoff modes of the fermion deter- 
minant is proposed. The level of suppression can be tuned by a new parameter /i 
that enters the theory. For fj, — » oo one recovers Wilson fermions whereas for /i — > 
one gets the quenched theory. For general p, the theory is similar to a Wilson theory 
with the effective cutoff ~ fi/a. 

By fixing /i at the appropriate value one can reduce considerably the UV fluctua- 
tions in the fermion determinant while preserving the interesting infrared modes of 
the theory. Thus one can clcminatc the need to separate and then devise algorithms 
for the UV modes. 

The fermion thoery described here can be easily extended to other input fermions 
such as staggered or chiral fermions. 



Appendix A: Proof of the result (1.7) 



I follow the analogous arguments as in Appendix C of [Hernandez et al, 1999|. 
Using the definition of the Wilson operator (1.1), a straightforward computation 
yields: 

A(U) = D{U)^ W D{U) W = 



+m 2 t - ma E M =i V M 

+|E^{-7M7.[V M ,V l ,]+a7 M [V M ,A i/ ]} 
Aa 2 [V^V^{x) = 

+ [U(x) /x U (x + ae^ v - U {x) v U (x + ae v )^]ip(x + ae M + ae v ) 
+3 similar terms 



where 



and 



2a 3 [V tl ,A„]ij(x) = 

+ [U(x)fiU(x + ae M )„ — U (x) v U (x + ae„) M ]^(x + ae M -I- ae u ) 



+3 similar terms 
Note also that: 



U(x)fj_U(x + ae M )„ - U(x) v U{x + ae v ). 



t-U v 
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and 

||a 2 [V M! VJ|| 2 
Therefore, I obtain: 



in 



|a 2 [V M ,A, 



2111- 



\\A(U)\\ 2 < \\D(t) w \\l 
One can easily show that: 



2^ 



|l-Up|| a 



||1 - U P \\ 2 < p - U P \\ F = v/6[l - P(U)] 



where ||.|| F is the Frobenius (Euclidean) norm of a matrix. From the definition (1.5) 
and ||D(1) W || 2 < m + f , \\SA(im 2 < 2 \\A(U)\\ 2 (assuming that \\A(U + SU)\\ 
||A(Z7)|| 2 ), I obtain the result (|L7|) with: 



< 



Cl 



2(m+ -) 2 , c 2 
a 



36V6 



Appendix B 

To prove the statement I use similar algebraic manipulations to those used else- 
where in a different context [Borici, 1999a|. Here they appear in greater detail. 
I let the lattice spacing in the fifth direction to be: 

L 5 

a 5 = — 
n 

where n is the number of lattice points in the fifth direction. The approximate 
fcrmion measure density can be defined by: 

det(l-T") 
Un{Hw) ~ det(l + T«) 
with T being a classical transfer matrix: 

T = 1 — a$Hw 

It is easy to see that for small lattice spacing 05, cu n (Hw) — > u(Hw)- It is only 
necessary to show that: 



det(D (5) 



w 



dct(l-T r ' 



where T ~ * T as as — > 0. The right hand side can be realized for example as the 
determinant of the following n x n block matrix: 

/ 1 -f \ 

11 '■• 



(4.1) 



V-T 



-T 
1/ 



("5") 

where the sign of the left lower corner reverses if the boundary conditions of D- 



change from periodic to antiperiodic. Therefore, the ratio of the two determinants 
will be given by: 



det(l-f") det(l-T n ) 
det(l + T") ^ det(l + T n ) 



12 



ARTAN BORigi 



(5) 

One must now calculate T from D w . The fermion matrix can be written as an 
n x n block partition in the fifth dimension: 

fa b D w -1 P+ P_ 
P_ a 5 D w -1 



D 



(5) 
w 



V 



P+ 

P- a 5 D w - tj 



where P± are the usual spin projector operators in the fifth direction. I multiply 
the above matrix from the left with the following permutation matrix: 



\ 



P. 
P+J 



and I obtain the following result: 

(a b P+H w -l a 5 P-H w + l 

a 5 P + Hw — 1 



75 



\a 5 P-H w + 1 



a 5 P„H w + 1 
a 5 P+H w - 1/ 



Comparing this matrix to that containing T (4.1 ) I arrive at the following expression 
for the transfer matrix: 

f= l )r-n-(t + a 5 P-H w ) 

1 — a 5 P + Hw 

which goes to T^ 1 for small lattice spacing 05. Since T and T -1 are equivalent 
transfer matrices the proof is concluded. □ 



Appendix C: Locality in the presence of gauge fields 
To show locality of D it is sufficient to show that the matrix- valued function /(.) 

, m £ N be defined by: 



(p.3) is a local operator. 

Let p— norm of the vector v £ 



NI P = (X>H 5 
fc=i 

with p > 0. The induced matrix norm of A £ R mxm is defined by: 

\\A\L = max \\Av\L 
p IMI P =i p 

If x, y are two lattice points in four dimensions, I will refer to the matrix element 
(x,y) of the matrix function / as f{x,y). Note that these elements are in fact 
matrices which carry spin and color indices, which I suppress for clarity. 

I will show in the following that: 

(4.2) \\f(x,y)\\ 2 < cxcl 11 ^" 1 , ci > 0, < c 2 < 1 

for any fi > and four dimensional lattice points x, y. 
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In order to prove (4.2) I need first to approximate the function /. From Appendix 
B one can infer that a five dimensional formulation gives a fermion measure density 
proportional to 1 — T n . Therefore one can write: 



D 



1 1 

-75- 



rpri 



a 1 + T n 

This suggests the following approximation to the function /: 



f(z 2 ) « f n (z 2 ) := 



— ) 

2n ' 



2n 



(1 - -^) 2n 
^ 2n> 



In' 1 V 2n> 

The right hand side can be expressed as a partial fraction by computing its poles 
and corresponding residues: 

1 ^ 1 



(4.3) 



Now I consider the matrix valued function /„ with the matrix aH ^ v substituted 
for the variable z and show first that it is local for any approximation order n. 
To simplify the notation, I call A k the matrix: 



(4.4) 

and write: 



A k = { 



,aH 



— — rcos 2 — (k ) + sin 2 — (k ) 

2np ' 2n 2 J 2n 2 J 



fn[( 



W .21 



2n 2 ^ A k 

fe=i K 



I use a geometric series to approximate the inverse of A k in order to keep the 
discussion simple (however note that a Chebyshev approximation is more accu- 
rate and is likely to decrease the exponent base C2 of (4.2) in the final result 
|Hcrnandcz et al, 1999fl ). I may write: 

(4.5) ' ' ' 



A k ||Ak|| a l-(l- 



11-4*11 



-) 



and then define: 



Pk = 



l- 



\M 2 



From the definition of the 2-norm I obtain: 

1 



Pk = 1 - 



\Ak\\ 2 A 



-ii 
fc I 



n{A k ) 



where K(A k ) is the condition number of A k . It is clear that p k , k = 1, . . . ,n are 
non-negative (k(A/s) > 1). 

Let o i < (72 be the extreme singular values of aDw- Then p k can be written as: 

& 2 ™ 2 £;(k-l) + sm 2 £(k-l) 



Pk = l - 



(JS_)2 cos 2 JL(k - I) + sin 2 J-(k 

V 2n/j. / 2n v 2 / 2n v 
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But since tanz > z for < z < ^ (the range of tan argument values for k 
1, . . . , n) I obtain: 



1 - 



Pk < 



[MTT(ft-i)] 2 



Clearly, I have: 
(4.6) 



p k < pi < 1, fc = 2, 



Now I may expand the right hand side of (4.5) in geometric series to obtain: 
1 

A~k 



Further, I let the matrix Mi. to be: 



' [i+(i- 1 r^ i -) + (i- Ak ^ 2 



fella 



L4 



fella 



Af fe = (1 - 



l^fella 



such that the matrix elements of /«(•) can be given by: 

1 " 1 



,l^4feho 

fe=l 11 " ll2 Z>0 



Since j4fc(a;,?/) vanishes for those lattice points x,y such that \\x — > 2a, one 
may conclude that: 

(4.7) Mj[(a:,tf)=0 1 for \\x-y\\ 1 > 2la, 1 = 1,2,... 
In general I have: 

(4.8) \\M l k (x,y)\\ 9 <1, 1 = 1,2,... 



Using (^7D and (|4.8|) , I obtain 



||/„(^ y )|| 2 <Pf l|a: - y|l Vn[(^) 2 ] 



ll/n^y)!!^^ 11 ^" 1 

^2 



But zf n (z 2 ) < 1 for z > 0. Hence finally: 
(4.9) 

Since the left hand side is uniformly bounded the result also holds in the limit 
n — > co. Hence the hypothesis (4.2) is proven with the constants c\, C2 given by: 

1 

J- rrr 



(4.10) 



Ci := — , and C2 := 

02 



1 _|_ M 2 ^ 2 



□ 

Remark. If the Wilson Dirac operator is singular (at o\ = 0, or at the "critical 
hopping parameter" ) the locality of the theory is guaranteed solely by the positivity 
of jl. 
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